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- - I. INTRODUCTION 

o ■ 

It is widely accepted that Einstein's general theory of relativity is an effective four-dimensional theory at lower 
energies and that it needs modifications at high energies. In the cosmological context, this break down of the theory 
results in a modification of the Fricdmann equation with important consequences in the dynamical evolution at 
early times (see for example, [H, Hj). There has been several proposals for cosmological scenarios mostly inspired by 
brane-world models or by the presence of dark energy. All of them lead to modifications to the total energy density 
dependence of the Friedmann equation. Also, there have been studies using dynamical system approach, bringing 
together several of such corrections to the Friedmann equation, in order to study the conditions which lead to attractor 
scaling solutions 0|. 

In the dark energy scenarios, the energy density of the scalar field, responsible for the dark energy, is im por tant 
VO \ onl y a t late times and can account for the acceleration of the universe 0, H, H, 0, H, H, OH EH, El, Eil, El, EH lid. E3, 
EH. Il9l. I20I l2ll I22I [23j . Therefore, scaling solutions (i.e. solutions for which the energy density of the scalar field and 
the perfect fluid scale in the same way), which appear as late time attractors of the system of the evolution equations 
can play an important role in modeling dark energy scenarios (24| ■ The existence and behavior of these models in the 
scaling regime depends on the proposed modifications to the Friedmann equation, as well as on the form of the scalar 
\ field potential. 

In this paper we analyze asymptotic evolution of modified cosmologies with a scalar field and a barotropic perfect 
fluid. The modification of the standard cosmology is parameterized by a function of the total energy density, in the 
line of Copeland et al. ||. The analysis of the asymptotic behavior of the equations is performed using the techniques 
^jjy of the renormalization group (RG), which has emerged as a powerful method for doing global and asymptotic analysis 
of ordinary and partial differential equations (Illinois group (25| , Bricmont et al [2(| , Caginalp [27| , Moise and Ziane 
[111 , for a reformulation of the RG approach in terms of the classical theory of envelops see Kunihiro, [2i| , the theory of 
perturbations of an isotropic universe with dynamically evolving Newton constant and Cosmologial constant [30l l3l| , 
and [H, HH for application of RG techniques to General Relativity). We shall show in this paper that the RG 
technique provides a very effective method for obtaining and analyzing such solutions. 

The plan of the paper is the following: In section II we describe the basic equations governing generalized cosmology. 
This is followed by a section where we give a very brief introduction to the RG method and apply this technique 
to the generalized cosmological models. Next section deals with a study of the exact scaling solutions, which are 
the fixed-points of the RG transformations. The stability of these exact solutions under linearized perturbations is 
analyzed in section V. We conclude with a brief summary. 
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II. GENERALIZED COSMOLOGY 



Extending their earlier analysis of FLRW models using dynamical system approach 0], Copeland et al. showed 
that the modified cosmological models, such as Randall-Sundrum, Shatnov-Shani, etc. can be recovered as scaling 
solutions by modifying the Friedmann equation 



H* = ^LpG 2 (p). (2.1) 



For G(p) = constant, we recover the standard FLRW cosmology. Here p is total energy density, and 777,4 is the 
four-dimensional Planck mass. 
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The model of universe considered here contains a barotropic fluid p 7 = (7 — l)p 7 , with adiabatic index < 7 < 2, 
and a scalar field with energy density 



1 



(2.2) 



where V{<j>) is the scalar field potential and p 7 
equations 



H = - 



4ttG 



-3H 



p. From the Einstein equations we get the following system of 



(2.3) 



G + 2pf ) [(l-^)^+l(p-V) 
(l--±)^ +1 (p-V) 



P = - 

p 7 = —3jHp 7 

The evolution of scalar potential V and the generalization function G depends on their respective functional depen- 
dence on </>, and p , 



\dp 

and which are for the moment left arbitrary. The generalization function G modifies the Hubble parameter, changing 
the cosmological expansion/contraction rate (H = a/a), and that correspondingly alters the evolution of all the 
physical variables in the system. 

Before applying the RG procedure to the above set of evolution equations some comments on scaling solutions are 
in order. It is well known that the scalar fields with exponential potential admit scaling solutions. The form of the 
potential changes in the case of generalized cosmologies; Copeland et al. in ref. [|[ gave a necessary and sufficient 
condition relating the scalar field potential and the generalization function G in order to get scaling solutions. Instead 
of repeating their procedure to get the required condition, we give below an alternative derivation to obtain an 
equivalent condition. 

Let's introduce a new variable 



Essentially we want to use the scale factor a(t), which is related to Hubble parameter, instead of the usual coordinate 
time t, as our evolution parameter. In the equation set (|2.3p . the equation giving evolution of the energy density of 
the fluid can be integrated trivially in the new parametrization as 



P7 — Po7^ 



(2.4) 



where p Q1 is an integration constant. The argument leading to the desired relation between V and G is now straight 
forward. Restricting ourselves to scaling solutions, i.e. both the scalar field density and the fluid density scale in the 
same way, we have 



P<p ~ Pi 



P<t> = Po<p3 



p = p a x 



(2.5) 



where p Q — p $ + p Q1 . Using the evolution equation of the total density along with the scaling condition above, the 
form of potential can now be determined as 



V = — - — Po^>X. 



(2.6) 



Also, using definition of , 



the kinetic energy term for the scalar field can be recovered as, 

<i> 2 = ipo<px. 



(2.7) 



On the other hand, we have 



dx 



(2.8) 
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Taking the square of the above equation, substituting (|2.ip and (|2.7p . and bearing in mind that now function G 
depends on x, we get 



TO - = ±/_£_. (2.9) 



3m 4 2 y p O J xG(x) ' 

Once the function G is fixed, the above equation gives <f> as a function of x. Now taking the inverse of this function 
and substituting it in (|2.6[) the form of potential corresponding to the specific choice of G is obtained. The condition 
given above is equivalent to the one obtained by Copeland et al. [3j. Therefore, condition for the scaling solution, 
relating generalization parameter G with potential V, can be easily recovered starting with the scaling ansatz, that 
is, both the scalar field component and the fluid component scale in the same way. 

III. RENORMALIZATION GROUP APPROACH 

The renormalization group (RG) approach developed by Wilson (see for example, Wilson and Kogut [34[) has 
enjoyed tremendous success, and is seen as a broad philosophy rather than just a mathematical technique. In the past 
decade this philosophy is also applied with considerable success to understanding some very basic aspects of non-linear 
differential equations. In particular, equivalence of RG and intermediate asymptotics was shown by Goldenfeld et. al. 
(see for example, [25[). They showed that the so called anomalous dimension in RG theory are actually the non-trivial 
exponents appearing in the intermediate asymptotics. 

The scaling solutions which we are concerned with are exact solutions, and which exhibit self-similarity at an 
asymptotic fixed-point. Therefore, RG seems to be a natural choice for our purpose. The form in which RG is used 
here was developed by Bricmont and Kupiainen [261 ]. The basic idea is as follows: We do a finite time integration 
followed by a rescaling of variables in the problem. This scaling transformation together with evolution equations gives 
us the RG equations. And, the fixed-points of these equations are actually the scale-invariant solutions to differential 
(evolution) equations which we are interested in. 

We will apply now the RG method to the system of equations governing the generalized cosmology, given in the 
previous section. Our interest here is in the early or late time regime, where solution asymptotes to the form 



lim u{t) = t x u*(l). (3.1) 



Here u(t), refers to any variable of the former system of equations, namely H, p,<p,V,G. The value "one" in the 
argument of lt*(l), signifies the initial value of the u, and \ is t ne scaling exponent to be determined later. It is 
convenient to choose t — 1 as initial time. The RG method gives us a systematic procedure to determine the exponent 
X, and it*(l), as we have illustrated below for the case of generalized cosmological models. Apart from the large time 
decay, this RG procedure is also applicable to the finite time decay or blow up of solutions, where self-similarity is 
exhibited at the asymptotic fixed-point. 

Let us consider following scale transformations 

t^Lt, H -> L e H(Lt) = H L , p — > L a p(Lt) = pl 

L b 4>{Lt) = 4> L , V^L C V = V L , G — > L d G = Gl (3.2) 

We use here a number L > 1, as a parameter for scale transformations, and the quantities subscripted with index "L" 
are scaled quantities. The RG transformation TZl is defined as a map from one initial data set to another [32j . 

TZ L u(x, 1) = ml(x, 1). 

These transformations have the semi-group property TZr,n — TZ^n-i o TZl- Since the scaled quantities satisfy the 
original system of equations, this fixes the following exponents 

e = 1, a = 2b = c, d= 1 - | . (3.3) 

Moreover, letting t = 1 and L = t, now we can express solution at arbitrary time in terms of initial values 

H(t) = t^H^l), P(t) = t- a PL (l), j>{t) - rH L (l), V(t) = t- a V L , G(t) = ti^Gt (3.4) 

where quantities sub-scripted with L are constants. The exponent a is the anomalous dimension: the field equations 
alone do not fix it and it should be determined by initial or boundary conditions. For instance, when one chooses a 
particular model by fixing the function G, parameter a gets determined. 
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Note from Eq. (|3 .4|) that the evolution of <j> (as well as <ft) depends on the value of parameter a. This behavior, 
as we will see later, has consequences on the nature of the potential in the scaling regime. In particular it is the 
parameter a, which distinguishes FRW standard cosmology from other "generalized" models. 

Comments are now in order about nature of the potential V. Though the functional form of the potential is left 
free, clearly in the scaling regime it gets fixed, Eq. (|3.4[) . Moreover, the functional form of G also is not arbitrary. As 
argued in the previous paragraph, depending on the value of parameter a, the scaling regime can be divided in two 
categories; a = 2, which gives the usual FLRW models and a ^ 2, the generalized cosmology. 
a = 2 : From Eq. (l3~4")l . we get 

V{4>) = V exp(-20/0 L ), G(p) = G . (3.5) 

Where Vq = Vl(1) and Go = Gl(1) are constants. Therefore, in the scaling regime we have a exponential potential 
and the generalization function G is a constant, i.e., FLRW cosmology. 

a ^ 2 : In this case G and V have a power law dependence on p and 4>, respectively, in the scaling regime 

V{4>) = Vo^l & , G(p) = Gopf? , (3.6) 

where Go and Vq have simple dependence on the L sub-scripted quantities. 

We draw a comparison now between the scaling solutions in cosmology, discussed in the previous section, and the 
scale invariant solutions which arise as asymptotic fixed-points of the RG analysis. First note that we are using as a 
variable the total energy density p, instead of the density of the perfect fluid p 1 . However, from Eq. (|2.3[) we have the 
following scaling relation for the perfect fluid energy density p 1 ~ t~ a . It is easy now to see that the energy density 
of the perfect fluid and the energy density of the scalar field scale in the same way 



k 2 + 2V L 



Pi 2 p 



constant (3-7) 



1L 



which is the usual definition of scaling solutions in scalar field cosmologies. Therefore, in RG method we recover all 
the scaling solutions of the system as a subset of scale invariant solutions, provided that both components of energy 
density are non-vanishing. 

The large L in our notation means late time. Therefore, applying the RG transformation repeatedly we can recover 
the long time behavior of the solution. Defining a auxiliary parameter r through L — exp(r), the RG equations are 



dH T . _ dH T 



dr dt 



= H L - 


\^G L ( 




4=1 


I rnj [ 



dG 



dpL 



G L + 2p L ^\ [l-l)4>l+l{p L -V L ) 



dpL 



dpi 



dr ^ L ' dt 



= a PL - {3H L [(l - |) 4>\ + j(p L - Vl)] } <=i 



d4>L _ a ■ d4>L_ 
dr ~ 2 <?i dt 



I .„ ■ dV L \ 



b L J t =i 



(3. 



dV L _ ^ | dV L 
dr dt 



\d(p L J j t=1 



dG L 
dr 



1 - - G 



dG, 



dt 



1 - - G 



(Pi'- 



Which is a set of algebraic equations since the quantities on the right hand side are all evaluated at t = 1. Note that 
the system of equations above has no dependence on L; this is expected since the system (12.3[) is scale invariant. 



5 



IV. SCALE INVARIANT SOLUTIONS 

The scale invariant solutions appear as fixed-points to the RG equations. The fixed-points are those points which 
are mapped onto themselves by the RG transformations for any L > 1 

IZlu* = u* , 

and which is equivalent to 

Therefore, the equations determining the fixed-points of our system are 

dm dp* deb* dV T * , N 

-7^ = 0, -7^=0, -p = 0, -t^ = 0. (4.1) 
dr dr dr dr 

As discussed earlier in this paper, if classified in terms of the potential function of the scalar field there are two 
classes of solutions corresponding to the exponential and the power-law potentials; equivalently the system can also 
be classified in terms of the generalization function G. Moreover, these two different classes also differ in their value of 
anomalous dimension "a" and which, therefore, can also serve as an alternative way of classifying these fixed-points. 

A. Exponential Potential (a = 2) 

To facilitate a comparison with earlier works in FLRW scalar field cosmology we shall use p 7 , as a variable in this 
subsection instead of the total energy density p. Moreover, we work with the notation used in [7[, defining k = 8irG. 
The set of RG equations simplifies to 

H L - <j — {1P 1 L + 4>l ) 

' t=i 



dr [ 1 

dp 7 L 



dr 



5L 



2 PlL -i 1 {p lL H L } t=1 (4.2) 



Ml 



dV L \ 



dr "' \ d(j) L J t=1 

Therefore, the fixed-points for this system are 



^-0, ^ = 0, ^=0. (4.3) 



dr dr dr 

The fixed-points should also satisfy the Friedmann equation 

H* 2 = ^{p* 1 + \4>* 2 + V*). (4.4) 

From the second equation in Eq. set (|4.3j) . we see that there are two different sets of fixed-points: those with p* = 
and those with p* ^ 0. Let's start with the fixed-points characterized by p* = 0. This case can be further subdivided 
into two sub-cases. In the first case, corresponding to V* — 0, we obtain 

dW I , 1 [2 

3' ^ = °' * = ± kV3' 



v * = ^ I — ) H* = -, P ; = 0, 0*=±- A /-. (4.5) 



This solution corresponds to a massless scalar field cosmology. 

In the second case, V* ^ 0, the nature of potential in the neighborhood of the fixed-point is given by (|3.5I) . 
Comparing the potential (|3.5p with the exponential potential used in Q, V = Vo exp(— nXcb), the fixed-point is given 
by 

^ 4=°> H * = i?> ^ = ^(^- 1 V (4 - 6) 
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It is a exponential potential scalar field cosmology. 

For the fixed-points with non-zero perfect fluid p* ^ the value of Hubble parameter is determined H* = 2/37. And 
we have following sub-cases due the nature of the potential. One case with V* — 0, and the solution for fixed-point is 



2 

37' f " r 3k 2 7 

Which is a perfect fluid only cosmology. The second case V* ^ 0, the fixed-point is 

7-3A 2 • J_ 2(2-7 ) 
3 7 ' r7 k 2 7 2 9A 2 ' 9 3kA' 9k 2 jX 



H * = — > Pi = ^,2 nx2 ■ ^ = — > ^* = n..2-.x2 ■ ( 48 ) 



for A 2 > 37. 

The fixed-points obtained are the same as in given in Q, where the scalar field potential is assumed to be of 
exponential form from the very beginning, and the system of equations was written in terms of normalized variables 
that make the phase space of the autonomous system bounded. It is interesting to note that imposing scale invariance 
we obtain the same results without using normalized variables. On the other side, it is reasonable to think that, what 
are called scaling solutions should be scale invariance solutions, since in that case the scalar field energy density and 
the energy density of the matter scale as a power of the scale factor. Assuming this we have proved that the only 
potential that gives scale invariance solutions is of exponential type in FLRW cosmology. 



B. Power- law potentials (a 7^ 2) 

In deriving the expression for potential in the scaling regime we have used the scaling behavior of <j), i.e., 



implying 20* = (2 — a)(f>* at the fixed-point. 

Let's note also that, as in the standard cosmology, there are two classes of fixed-points depending on whether p 7 
vanishes or not. The quantity p 1 , however, is not a variable in our system of equations (|3.8[) (we have used p instead), 
nevertheless /j 7 satisfies the following equation 

= P-r( a ^ 3 7#)> 

which means that if fluid density is non- vanishing at the fixed point p* ^ 0, the Hubble parameter takes the critical 
value H* =0/37. 

To analyze the fixed-points when a ^ 2, we first consider V = 0. In this case the RG equations simplify considerably 

87rG 2 



giving the following set of fixed-points 



2ap* = 3ir[(2-7)(/>* 2 + 27p*] (4.9) 
a<j)* = 6H*<j}* 



{J ^Y. r-± (4.10) 



V247T7 2 G 2 1 y' ' 3 7 



( a 2 



The fixed-point (|4.10p corresponds to vanishing of the scalar field and (|4.11|) to the vanishing of the perfect fluid. 
Therefore, only one component survives in the scaling regime. Note that (|4. 1 1[) is not a limiting case of (|4.10p for 
7-2. 
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The fixed-points in the more general case with non-zero potential are governed by the following set of algebraic 
equations 

H* = —^-p*° 
6m\ 



8ap* = 3H*[(2 - 7 )(2 - a) 2 cj)* 2 + 8 7 (p* - V a (f>*^)} 



(4.12) 



H*6* = 



4aVn 



3{a-2) 2 



These cannot be solved in general however various special cases corresponding to different values of parameter 2) 
can be analyzed. 

First we consider the simplest scenario where potential V vanishes in the asymptotic regime due to vanishing of the 
scalar field itself (<p* = 0). This is unlike the previous case where potential was identically zero (due to, say, Vq = 0) 
to begin with. We recover following fixed-point for vanishing of the scalar field 0* = 0, 



H 



= — , P = — 



37' 



V247T 7 2 G2 



, V* 







(4.13) 



This is similar to the previous case, except the fact that the range of anomalous dimension a is restricted to \a\ > 2, 
and that the cosmology is driven towards vanishing potential which is non-zero to begin with. 

The fixed-points corresponding to a non-vanishing potential in the scaling regime can be divided in two families. 
The one corresponding to vanishing fluid component in scaling regime is given by 



H* = ±* 



'8^4 
3m 2 P ■ 



P = 



3tt(2 - a) 4 G 2 
2a 2 m 2 



2(a-l) 



V* = V cj)* — . 



(4.14) 



The plus and the minus signs correspond to expanding and contracting models, respectively, and the <t>* in the equation 
above is given by the solution to following algebraic equation (o^ 1,2) ; 



±1 



/247rG 2 



a 2 ml 



The case a = 1 needs to be treated separately and the fixed-point is of the form 
H 



( 2 -a? n 
+ — = 0. 



8F 



(4.15) 



3ml 



/37tg; 



I ( m 4 



+ Vo), 



,*2 



2ml 



V* = V <j)*~ . 



2m 2 4 V 127rG o 

The other fix-point corresponding to non- vanishing of both the scalar field and the perfect fluid is of the form 

(2- 7 )(2-a) 



H* = — , 

37 



47 Vb 



( cp-m 2 



247T7 2 G^ 



V* = Vol 



(4.16) 



(4.17) 



For a — 2 system reduces to usual FLRW models with exponential potentials which were analyzed in the previous 
subsection. 



C. Stability 

The scaling solutions studied in the previous section arise as fixed-points of the RG equations. The stability in 
time evolution around a self-similar solution is stability of RG flow near the fixed-point. We do a linear perturbation 
analysis here, considering perturbations in the initial data space staying close to the fixed-points and searching for 
relevant modes, i.e., modes which make the flow diverge from the fixed-points. 

Let us define following perturbation quantities 



H = H* + SH, p = p* + 5/?, 4> = 4>* + 5<p, 4> = <p*+ 6<j>, V = V* + 6V, G = G* + SG 
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Where all perturbed quantities are small (<C 1)- The Friedmann equation can be linearized as 

4-7T 



H*5H 



3m| 



G*[G*Sp + 2p*SG] 



(4.18) 



and can be used to relate variation in H with other variables. However, this reduction in variables can make the 
equation set much more complicated and for the moment we will work with the full set of variables. 
The linearized RG equations take the following form 



dSp 
~dT 



= (a-3jH*)Sp + 3{j-2)H*(, 



3jH*6V - 



3 r 



(2- 7 )(0*) 2 + 2 7 (p* 



V*) 



SH 



dr 



- 3iT 



3^* SH 



2a V* 
(a - 2) (0*)2 



2a 



(2 -a) 



-SV 



dr 
dSV 
~~dT 

dSG 

~~dT 



(1 



2a V* 
(a -2)~ft> 
(a -2) 



2a V*<j>* 



(a - 2) cj)* 



SV 



2a 

+2jV*} dp 



3H* 



(a-37ff*) + ~ ( 7 "2)(0*) 2 + 2 7 y* 

Z p L 



<5G + - 



3 G*iT 



7G*-^ [( 7 -2)(0*) 2 + 2 7 y* 
2 p* L 



«? + 3(2- 7 ) 



2 

H*G 



(7-2)(<n 2 



We now find the normal modes with following ansatz 

= /exp(wr) 



(4.19) 



where / is a representative variable for H, p, <fi, </>, V and G, and over bar signifies quantity to be constant. With this 
ansatz and also simplifying system using (|4.18p we can now look for any relevant modes. 

Of all the fixed-points listed in the previous section the one given by (|4.17[) is of particular interest since this 
corresponds to late time regime where both scalar field as well as perfect fluid component are non-vanishing. Taking 
these particular fixed-point values the stability of the system depends of the following matrix 



i. 



q(7-2)(a-2) ,* 

27 * 



6x7(2 — a) q*2i~ 



a(2-7) 
2 7 



q(q-2)(7-2) 
47 



1^27 



12 ^ {2 .- a) G*p*(h* 

rt.mr. " " 



UJ + 1 



q(q-2)( 7 -2) ^ 
47 



(2- 7 )(q-2)* G*4>" 
47 p* 



q(o-2)-(7-2) ,» 
82 







127T7^(q-2) ^2 „* 



V 



6^(2-0)^*3 

9 2 '- ' 







(q-2) g* 
2 P * 



w ■ 



G* V / 



The values of arc determined by vanishing of the determinant of the above matrix. The two eigenvalues are zero 
and the remaining three are determined by solution of a cubic equation 



(l-a+- 

7 



247ry 2 , 

2 2 G P 



0(2-7) 67r( 7 -2)(a-2) 2 c , 2 2 



24^7(7 + a — 7a) 



-.2 ^-v-, 2 



7 

6tt(7 - 2) 



am: 



(a-2) 2 0* 2 +4 7 p* 



. 



(4.20) 



Solutions of this equation and therefore criterion for stability of the fixed-points depends on the value of the anomalous 
dimension a and and other constants in the problem. 



V. SUMMARY 



We have analyzed in this paper long time behavior of the cosmological equations By using the Renormalization 
Group method, paying special attention to the scaling solutions in two different cases. The first case of standard 
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FRLW cosmology has been widely studied and it serves us to illustrate the method we use. Although the system in 
this simple case can be analyzed by means of dynamical system technique, this analysis depends upon the choice of 
normalized variables which are not easy to find. The RG method avoids this problem: the scale invariance solutions 
found describe all the scaling solutions. 

We also applied RG method to the generalized cosmologies, which imply modifications of the Friedmann equation. 
The RG equations give the scaling solutions of the system and we have shown that in the scaling regime the potential 
of the scalar field and the function which describes the modifications to the Friedmann equation have a power law 
dependence on their respective variables. 

The RG method is a simple but powerful tool to investigate scaling solutions for extended cosmologies. We have 
proved that to get scaling solutions, during the dynamical evolution, the scalar field potential should be either power 
law or exponential and the function G, that generalizes Friedmann equation, should be a power law of the total 
density. Only the fixed point given by Eq. (|4.17p represents a scaling solution, whereas the rest of fixed points give 
cosmologies with either vanishing fluid or vanishing scalar field. The stability of the scaling solution fixed point is given 
by Eq. (|4.20[) , and it depends on the particular model under study through the constants appearing in that equation. 
However, the most relevant of these constants is the anomalous dimension "a". Finally, it would be interesting to 
study the evolution of a larger variety of field models as K-essence, phantom, quintaessence, etc. (see [35| for a review) 
with this technique. 
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